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The physics of turbulence remains one of the most challenging problems in fluid dynam- 
ics. Although more than a century of effort has been devoted to it, a lot of fundamental 
issues axe still unresolved. This is particularly so in the case of turbulence in high speed 
flows because of the increased number of possible mode interactions due to compressibility 
effects. For example, the cubic non-linearities in the momentum equations allow the vor- 
ticity, acoustic and entropy modes to interact with each other. The dynamics are further 
complicated by the possible existence of non-stationary shocks and/or eddy shocklets. 

In this paper, several direct simulation^ of 3-D homogeneous, compressible turbulence 
are presented with emphasis on the differences with incompressible turbulent simulations. 
A fully spectral collocation algorithm, periodic in all directions coupled with a 3rd order 
Runga-Kutta time discretization scheme is sufficient to produce well-resolved flows at Tay- 
lor Reynolds numbers below 40 on grids of 128x128x128. A Helmholtz decomposition of 
velocity is useful to differentiate between the purely compressible effects and those effects 
solely due to vorticity production. In the context of homogeneous flows, this decompo- 
sition is unique. Time-dependent energy and dissipation spectra of the compressible and 
solenoidal velocity components indicate the presence of localized small scale structures. 
These structures are strongly a function of the initial conditions. We concentrate on a 
regime characterized by very small fluctuating Mach numbers Ma (on the order of 0.03} 
and density and temperature fluctuations much greater than Ma 2 , This leads to a state in 
which more than 70% of the kinetic energy is contained in the so-called compressible com- 
ponent of the velocity. Furthermore, these conditions lead to the formation of curved weak 
shocks (or shocklets) which travel at approximately the sound speed across the physical 
domain. Various terms in the vorticity and divergence of velocity production equations are 
plotted versus time to gain some understanding of how small scales are actually formed. 
Possible links with Burger turbulence are examined. 

To visualize better the dynamics of the flow, new graphic visualization techniques have 
been developed. The 3-D structure of the shocks are visualized with the help of volume 
rendering algorithms developed in house. A combination of stereographic projection and 
animation greatly increase the number of visual cues necessary to properly interpret the 
complex flow. The presence or absence of shocks is automatically detected by monitoring 
of the minimum and maximum divergence of the velocity field over the physical domain. 
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NAVIER-STOKES EQUATIONS 
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PHYSICS 


g 

a 



H-J 

S3 

.1 

tu 

b- 

$ 

a 


C/5 

o 

bO 

0) 

05 

£ 

05 

T5 

li 


■s 

o 

CO 

o 

o 

M 

Pn 

"2 

o 

a 

• 

• 

• 

• 

• 


0 ) 

a 

05 


154 


INITIAL CONDITIONS (1) 
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INITIAL CONDITIONS (2) 
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NUMERICAL METHOD 
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turb3d/rur»1 32 
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K.E> dissipation 


turb3d/ run 1 32 
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OIUERGENCE OF UELOCITV 

turb3d/run13Z, tsS.3, it=600 
Re=150, Ha= . 028, k0=1.85, urns:.1, chis.07, rho,T=10* 
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CONCLUSIONS 

• W eak shocks can be generated by an initially random 2-D velocity and ther- 
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